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We address the question of how stock prices respond to changes in demand. We quantify the 
relations between price change G over a time interval At and two different measures of demand 
fluctuations: (a) $, defined as the difference between the number of buyer-initiated and seller- 
initiated trades, and (b) Q, defined as the difference in number of shares traded in buyer and seller 
initiated trades. We find that the conditional expectations {G)n and {G)<s> of price change for a 
given !] or $ are both concave. We find that large price fluctuations occur when demand is very 
small — a fact which is reminiscent of large fluctuations that occur at critical points in spin systems, 
where the divergent nature of the response function leads to large fluctuations. 
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Stock prices respond to fluctuations in demand, just as 
the magnetization of an interacting spin system responds 
to fluctuations in the magnetic field. Periods with large 
number of market participants buying the stock imply 
mainly positive changes in price, analogous to a mag- 
netic field causing spins in a magnet to align. Here, we 
quantify how price changes depend on demand fluctua- 
tions [|J-g|, and find a strikingly non-linear relationship 
with a specific functional form which is not altogether un- 
like the dependence of magnetization on field strength. 

Fluctuations in demand arise from buy and sell orders 
of market participants. To quantify fluctuations in de- 
mand, we distinguish buyer-initiated and seller-initiated 
trades defined by which of the two participants in the 
trade, the buyer or the seller, is more eager to trade. 
When such a distinction does not exist, we label the trade 
as indeterminate. We identify buyer and seller initiated 
trades using the bid and ask quotes Sb(*) and iSa(£) at 
which a market maker is willing to buy or sell respec- 
tively @). Using the mid-value S M {t) = (S A (t)+S B (t))/2 
of the prevailing quote (^-0] , we label a trade buyer initi- 
ated if S(t) > Su{t), and seller initiated if S(t) < <Sm(£)- 
For trades occurring exactly at S-M.{t), we use the sign of 
the change in price from the previous trade to determine 
whether the trade is buyer or seller initiated, while if the 
previous trade is at the current trade price, the trade is 
labelled indeterminate ||||]. Accordingly, for each trade 
i, we define the variable 



1 (buyer initiated) 
(indeterminate) 
— 1 (seller initiated) 



(1) 



Wc quantify demand fluctuations by analyzing two 
quantities: (a) the number imbalance (difference be- 
tween the number of buyer-initiated and seller-initiated 
trades P,fol in a time interval [t, t + At]), 



N 



(2a) 



and (b) the volume imbalance (difference between the 
number of shares traded in buyer-initiated and seller- 
initiated trades in a time interval At), 



n = ^At(t) 



N 

E 



(2b) 



where qi is the number of shares traded in trade i, and 
N = N& t (t) is the number of trades in At. 

To choose a time scale in which to analyze the de- 
pendence of price fluctuations on demand, we first com- 
pute the correlation functions [Fig. |l| ($(t)G(t + r)) and 
(il(t)G(t + t)) and find significant dependence at r = 0. 
For |t| > 0, both correlation functions decay rapidly, and 
cease to be statistically significant beyond r w 15 min — 
thereby setting a short time scale for the response of price 
changes to fluctuations in demand. 

Next, we shall examine the relationships 



(G) n = E(G|fi) 



(3a) 



(3b) 



which give the equal-time expectation value of price 
change G(t) for a given <I?(i) or il(t). Figures |^(a) and (b) 
show (G)$ and (G)n for 5 typical stocks for At = 15 min. 
We find that both (G)$ and (G)q are nonlinear, display- 
ing concave curvature with increasing 4> and fl [pd-p^[, 
and 'flattening' at large values [Q. 

Figure ||(c) shows the average behavior of {G)$ for all 
stocks. We find that (G)$ is consistent with the func- 
tional form 



(G)$ = A tanh(^i$) . 



(4) 



where Aq is a constant that denotes the level of 'satura- 
tion', and Ai determines the average price change for unit 



1 



change in <E>. In the case of a spin system, the saturation 
at large values for the analogous curve — magnetization 
vs. field — is due to the fixed number of spins. The 
apparent saturation of (G)$ is surprising in the present 
context, since there is no clear upper limit either on the 
price change, or on the number of trades. We find that 
(G}$ for a range of At, also displays good agreement with 
Eq. (@) [Fig. |(d)]. 

We next focus on (G) n [Fig. |(e)]. We find 
that the function (G)n, like (G)$, is consistent with 
Eq. (|) However, near ft = 0, (G) n shows 

not a strict linear behavior for small Q as we expect for 
tanhfi, but rather a power-law (G)n ~ [Fig. ||(e)]. 
We find that 1/6 depends on At [Fig. |(f)]: 5 w 3 for 
At = 5 min and 5 ~ 3/2 for 15 min, and 5 — > 1 for larger 
At, agreeing well with tanhSl 

Next, we analyze the dependence of the number of 
trades N on demand fluctuations to quantify how large 
volume imbalances generate trades. Figure ||(a) shows 
that the equal-time expectation value (iV)$ shows a lin- 
ear increase with The dependence of N on volume im- 
balance fl is nonlinear; (N)q displays a "cusp" at CI = 
followed by a sharp increase and saturation at large val- 
ues [Fig. |(b)] . 

In spin systems, the amplitude of spin fluctuations is 
related to the susceptibility, which quantifies the response 
of the system to fluctuations in the magnetic field. In our 
problem, a certain change A$ in demand $ (analog of 



the field) causes a response 



8<f> 



A$, which we find to 



be largest at $ = [Fig. |2j, suggesting that the non- 
linear shape of (G)$ can give rise to large fluctuations 
(large "volatility" [0) when $ is small. The average 
amplitude of fluctuations in G = 53i=i $Pi is given by 
the variance 



x 2 = (6p*) - {5 Pi f 



where 8pi is the price change due to trade i and (...) 
denotes the average computed over the interval At. Fig- 
ure ||(a) shows that (%)$ displays large values near $ = 
and a rapid decay for increasing $. Figure [|(b) shows as 
a function of <E> the number of events with price change 
|G| > 5 standard deviations. Interestingly, we find that 
a majority of the large events occur at $ = 0, consistent 
with previous empirical results [ 
power-law distribution of price c 

from x- Our findings are reminiscent of phase transitions 
in spin systems, where the divergent behavior of the re- 
sponse function at the critical point (zero magnetic field) 
leads to large fluctuations [ pp[ . 
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FIG. 1. Cross correla- 

tion functions (($(t)G(t+r))-($(f))(G(f)))/a G (T* (open cir- 
cles)and {{fl(t)G(t + T)) - (fl(t)}(G(t)))/crG<m (closed circles) 
computed using 5 min time series for <3>, fl, and G. We find 
short-range time dependence which after w 15 min reaches 
noise levels (dashed lines). 
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FIG. 2. (a) Conditional expectation {G)<s> of the price 
change for a given value of $ for 5 typical stocks over a time 
interval At = 15 min. Both G and <3? are normalized to have 
zero mean and unit variance, (b) Conditional expectation 
(G)n for the same 5 stocks as in part (a). We normalize G to 
have zero mean and unit variance. Since il has a tail expo- 
nent £ = 3/2 which implies divergent variance, we normalize 
O by the first moment (|fi — (0}|). (c) (G)$ averaged over all 
116 stocks studied. The solid curve shows a fit to the function 
J 4 tanh(y4i $), with A = 0.88 ± 0.01 and Ai = 0.38 ± 0.01, 
where the fit is performed with tolerance = 0.01. (d) Same 
as (c), on a log-log plot for $ > (filled symbols) and $ < 
(empty symbols) for At = 15 min and 195 min (shifted verti- 
cally for clarity). The solid curves show fits to Ao tanh(Ai $), 
which agrees well with the data, (e) Conditional expectation 
(G)n averaged over all 116 stocks. We calculate G and il 
for At = 15 min. The solid line shows a fit to the function 
Bq tanh(_Bif2). (f) (G)n on a log-log plot for different At. For 
small «, (G)n ~ tt 1/s . For At = 15 min find a mean value 
1/5 = 0.66±0.02 by fitting (G)n for all 116 stocks individually. 
The same procedure yields 1/5 — 0.34 ± 0.03 at At = 5 min 
(interestingly close to the value of the analogous critical expo- 
nent in mean field theory) . The solid curve shows a fit to the 
function Bo tanh(_Bi 57). For small Q, Bq tanh(Bi f2) ~ Q, 
and therefore disagrees with (G)n, whereas for large SI the fit 
shows good agreement. For At = 195 min (~ day) (squares), 
the hyperbolic tangent function shows good agreement. 
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FIG. 3. (a) Conditional expectation (iV)$ of the number of 
trades for a given $ averaged over all 116 stocks, shows ap- 
proximately linear behavior with increasing $. (b) {N)n av- 
eraged over all 116 stocks shows strikingly nonlinear behavior. 
The solid line shows a fit to the function Go — C\ exp(— C2O,) 
(which has the same large fl behavior as a hyperbolic tan- 
gent). For both parts (a) and (b) We calculate G, $ and Q 
over At = 15 min. Both $ and G are transformed to have 
zero mean and unit variance, whereas 51 is normalized by its 
first moment. 
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FIG. 4. (a) Conditional expectation (x)*> where x is calcu- 
lated using Eq. (|^), shows large values near <& = and decay 
for increasing The solid lines show a fit to the function 
Dosech 2 (Di$). (b) Number of events with |G| > 5 standard 
deviations for a given $ shows large values at $ = 0. 
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